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ABSTRACT 

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to ex- 
plore the specific star-formation activity of galaxies and their evolution near the peak 
of the cosmic far-infrared background at 70 and 160|J.m. We use a stacking analy- 
sis to determine the mean far-infrared properties of well defined subsets of galaxies 
at flux levels well below the far-infrared catalogue detection limits of SWIRE and 
other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to 
the far-infrared background at 70|_im and 160|^m. These long wavelengths provide a 
good constraint on the bolometric obscured emission. The large area provides good 
constraints at low z and in finer redshift bins than previous work. At all redshifts 
we find that the specific far-infrared luminosity decreases with increasing mass, fol- 
lowing a trend L Fm /M* oc Af f with = -0.38 ± 0.14. This is a more continuous 
change than expected from the iDe Lucia fc Blaizotl (|2007l ) semi-analytic model sug- 
gesting modifications to the feedback prescriptions. We see an increase in the specific 
far-infrared luminosity by about a factor of ~ 100 from < z < 2 and find that the 
specific far infrared luminosity evolves as (1 + z) a with a — 4.4 ± 0.3 for galaxies with 
10.5 < log 10 M„/Mq ^ 12. This is considerably steeper than the IDe Lucia &: Blaizotl 
(|2007| ) semi-analytic model (a ~ 2.5). When separating galaxies into early and late 
types on the basis of the optical/IR spectral energy distributions we find that the 
decrease in specific far-infrared luminosity with stellar mass is stronger in early type 
galaxies (j3 ~ —0.46), while late type galaxies exhibit a flatter trend (/3 ~ —0.15). 
The evolution is strong for both classes but stronger for the early type galaxies. The 
early types show a trend of decreasing strength of evolution as we move from lower 
to higher masses while the evolution of the late type galaxies has little dependence 
on stellar mass. We suggest that in late-type galaxies we are seeing a consistently 
declining specific star-formation rate a = 3.36 ± 0.16 through a common phenomenon 
e.g. exhaustion of gas supply i.e. not systematically dependent on the local properties 
of the galaxy. 
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1 INTRODUCTION 

A fundamental goal of modern astronomy is to understand 
the processes driving the formation and evolution of galax- 
ies. A key issue is the relationship between the assembly of 
galaxies, and the formation history of the stars within those 
galaxies. A galaxy can increase its stellar mass through the 
accrual of stars in a 'dry' merger where the merger does not 
trigger new star formation or directly through star formation 
triggered by a merger or some other process. The star forma- 
tion rates are determined by a variety of factors including the 
triggering mechanisms, the supply of gas and the feedback 
processes. The contribution from all of these three processes 
to stellar mass buildup in galaxies is subtle, and has been 
studied in numerous optical/near-IR photometric and spec- 
troscopic surveys. These surveys have demonstrated that, as 
we increase in redshift, there is a strong dependency on at 
least two parameters; galaxy mass, and local environment. 

Gala xy mass is thought to play a important role, at least 
at z < 1 |Cassata et alj|2007l ). At low redshifts (z < 0.2), 
ongoing star formation in massive galaxies is almost entirely 
absent. Extreme levels of star formation are found rarely and 
in many cases are triggered by interactions and mergers. 
Moderate star formation is probably triggered by internal 
processes ijOwers et al.ll2007l : iMelbourne et alj|2008l ). As we 
move to higher redshifts however, this picture changes. Over- 
all, the galaxy stellar mass function at high masses evolves 
fairly slowly up to z ~ 0.9, and then more rapidly up to 
at least z ~ 2.5, suggesting that th e majority of stellar 
mass assembly too k place at z > 1 jPozzetti et al.l 120071 : 
iFeulner et all 120071 ). Massive galaxies show little evidence 
for stellar mass assembly via either star formation or dry 
mergers at z < 0.7, while lower mass systems harbour on- 
going star formation at all redshifts, lending support to the 
idea of 'downsizing', in which more massive galaxies form 
most of their stars at hig h redshifts. Ther e is also some ev- 
idence for 'dry' mergers ()Bell et al.ll2006l ). Finally there is 
evidence that higher mass galaxies have lower specific star 
formation rates (sSFR, i.e the star formation rate per unit 
stellar mass) than do lower mass systems over a very wide 
redshift range, possibly up to z ~ 4, sug gesting that lower 
mass galaxies form stars more efficient ly ( Bauer et alJl200"S 
IFeulner et ai1l2005l ; Izheng et al ] |2007al ). although the sSFRs 
of massive galaxies appear to increase rapidly with increas- 
ing redshift. 

Local environment also has a significant effect. In the 
local Universe we see a distinct environmental segregation, 
in which galaxies in rich environments show much lower star 
form a tion rates than do ga l axies in the field (jCassata et al.l 
120071 : IZauderer et"al] l2007t ICooper et all l2007t ). At higher 
redshifts however this trend is reversed; at z ~ 1 the average 
star formation rate increases with increasing environmen- 
tal d ensity, as does specific star formation rate (jElbaz et al.l 
120071 ). and morphological evolution appe ars to be more rapi d 
in dense environments than in the field (jCapak et al ] l2007t ). 
Star forming galaxies at z ~ 1 are in richer regions than seen 
in the local Universe ()Farrah et al j200&lMagliocchetti et al.l 
120071 ). and there is evidence that star formation gradu- 
ally shifts to lower den sity regions from z ~ 1.5 to z=0 
|de La Torre et al"1l2007t ). thus providing a natural explana- 
tion for the local environmental segregation. Since environ- 
mental effects are so important we can only understand the 



universal properties when averaging over a representative 
sample of environments. 

An important consideration is the feedback from su- 
pernovae or AGN which is required to suppress star for- 
mation in semi-analytic models (SAMs). A growing consen- 
sus is that the models require AGN fe edback to suppress 
the star formation in massive objects jCroton et al.l [20061 : 
iBower et aI1l200cl) and this is in fact the dominant mecha- 



De Lucia fc Blaizotl l |2007f ). 



The specific star-formation rate is a particularly use- 
ful probe of these processes. It measures the ratio between 
the curre nt star-formation rate and the historically ave raged 
rate (e.g. iBrinchmann et alj|2004 IWalcher et alj|2008h . En- 
hanced star-formation rates arising simply from an increase 
in gas reservoirs achieved through 'dry mergers' will not af- 
fect the sSFR, which is mainly sensitive to the triggering, fu- 
elling and feedback i.e. the overall star-formation efficiency. 
Thus the sSFR provides some de-coupling of the merging 
and other phenomena. 

There remains however a significant problem with these 
studies. The discovery of a strong cosmic infrared back- 
groun d (CIRB) by COBE (|Puget et all 1 19961 : iHaiiser et all 
1 1998f ) . and its subsequent (partial) resolution into a huge 
opulation of obscured star-forming galaxies at z > 1 by ISO 
Rowan-Robinson et al. [ I l997l : iDole et al.ll200ll:lMann et al.1 



120021 : IVerma et al.1 120051). Spitzer dLe Floc'h et al 



and in the sub-mm (iHughes et al.l Il998l; Eales et al 



IScott et all 120021 : iBorvs et al.l 120031 : iMortier et al 



2005), 



.20001 : 
20051 ). 



demonstrated clearly that a large fraction of the total star 
formation at high redshifts is heavily shrouded in dust, and 
therefore impossible to detect at optical or near-IR wave- 
lengths. This applies to even mode rately luminous systems 
at z ~ 1 (e.g. iLe Floc'h et aLlliooE their figure 14). There- 
fore, optical/near-IR surveys to probe stellar mass assembly 
at high redshift miss a significant fraction of ongoing star 
formation. 

The Spitzer space telescope ijWerner et al.ll2004l ) has al- 
lowed us to make great progress in understanding these ob- 
scured systems. However, even the recent surveys that com- 
bine optical data with 24/xm data do not solve this problem; 
the peak rest-frame emission from most IR-luminous galax- 
ies is in the range 40-100/im, but at 1 < z < 2 the 24/im 
Spitzer band probes rest frame 8-12^im, a region which can 
be highly contaminated with a range of features (i.e. PAH 
emission, Si absorption etc). Spitzer does however have two 
longer wavelength channels, at 70 and 160/im. The resolu- 
tion and sensitivity of these channels are insufficient to di- 
rectly resolve the far-infrared (FIR) b ackground into t he in- 
dividual galaxies that produce it (e.g. lDole et al.ir2004l ). but 
techniques exist to alleviate this. One of these techniques is 
known as 'stacking'. 

To understand galaxy evolution it is thus essential 
to understand the obscured specific star-formation rate of 
galaxies (and massive galaxies in particular) over a represen- 
tative range of galaxy environments. In this paper, we use 
a stacking technique to measure the specific far infrared lu- 
minosity at 70 and 160um as a function of galaxy mass and 
redshift. To probe the highest stellar mass objects, which 
are rare, we use a large survey area. Although we don't ex- 
plore the variation with environment in this paper this large 
survey area also means we can be confident that we have 
covered a representative range of environments. 
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In Section 2 we discuss the samples. In Section 3 we 
describe our stacking technique. In Section 4 we present the 
results, which we discuss in Section 5, before concluding. We 
assume a spatially flat cosmology with Ho = 100 h km s _1 
Mpc - , h — 0.7, fi = 1, and fi m = 0.3 and all magnitudes 
are Vega magnitudes. 



(|Afonso Luis et al"1l2008l ) and find a good agreement (within 
the absolute calibration uncertainties). 

We note that the absolute SSC calibratio n is good to an 
accuracy of 7% for 70 n-m and 12% for 160um (|Gordon et al.l 
l200d IStansberrv et alJkood ). 



2 CATALOGUES AND MAPS 
2.1 SWIRE MIPS Maps 

The Spitzer Extra g alactic Legacy Sur vey (SWIRE, 
lLonsdale et al.l 120031 . I Lonsdale" et all I2004T I observed 49 
square degrees in six fields (CDFS, Elais-Nl, Elais-N2, 
Elais-Sl, Lockman and XMM-LSS) using the seven primary 
Spitzer imaging bands (3.6 to 160um). The SWIRE 70 urn 
and 160um maps used here w ere observed using MIPS 
medium scans as described in IShupe et al. | l|2008h . The 
70um maps were made from Basic Calibrated Data (BCD) 
images produced by the SSC pipelines, after applying time 
filtering and colum n filtering following the prescription in 
iFraver et al.l (|2006l ). and are part of our Data Release 4. The 
filtered BCDs were mosaic ked using the MOPEX package 
l|Makovoz fc Marleaull2003 ). The 160um maps were made 
by mosaicking the filtered BCD images produced by the 
SSC pipelines and correcting for a systematic 5" pointing 
offset. The maps are calibrated in units of surface brightness 
(MJy/sr). To bring the calibration in-line with the facility 
calibration appropriate for the "S13" pipeline processing 
(May 2006) we scale the maps up by a factor of 1.107 at 
70um. The 160um maps for CDFS, XMM-LSS and Elais- 
N2 were multiplied by 1.064 while those of Lockman and 
Elais-Nl were unchanged (having already been reprocessed 
with the latest calibration). I.e. our map calibration is 702 
MJy /sr per MIPS-JO un i t and 44.7 MJy/sr per MIPS-160 
unit (|Gordon et alj|2006l ; l&tansberrv et alj|2006l ). 

We estimate point-source fluxes, / (in mjy), by fitting 
the point source response function (PRF) to the map in- 
tensity / (in MJy sr _1 ). Our PRF is based on the default 
MIPS PRFs as recommended for use with Apex from the 
Spitzer Science Centr^E These PRFs have the same pixel 
size as the SWIRE maps i.e. 4" at 70um and 8" at 160um 
and are over-sampled by a factor of 4 i.e. l"(for 70um) and 
2" (for 160um). These PRFs are not identical to the ones 
used for the SWIRE catalogue extraction but the resulting 
difference calibration is less than 5% and smaller than the 
absolute calibration uncertainty. 

We define an effective beam size fibcam = f/Io 
| PdQ/Po (where P is the PRF, dQ, a solid angle and 
the zero subscript indicates the peak). We obtain fibcam = 
13.1 nSr at 70/mi and Slboam = 61.5 nSr at 160/^m. 

Finally, for consistency with the SWIRE catalogues we 
apply a colour correction to convert from the standard cal- 
ibration for a 10,000 K blackbody to constant uf v more 
appropriate for galaxies. We therefore multiply the 70 and 
160um fluxes by 1.09 and 1.043 respectively. 

We have checked our calibration by comparing the 
map intensity at the position of catalogued MIPS sources 



2.2 SWIRE optical/IR band-merged catalogues 

The fields Elais-Nl, Elais-N2, Lockman and CDFS currently 
have the most homogenous and well understood SWIRE 
data. The catalogues we are using are those that were re - 
leased as part of SWIRE Data Release 4 (|Surace et alj2004h . 
All IRAC and MIPS catalogues have been "bandmerged" i.e. 
independent catalogues from the seven different bands (3.6, 
4.5, 5.8, 8.0, 24, 70 and 160um) have been cross-matched 
to produce one master catalogue. In Elais-Nl and Elais-N2 
we have 5-ban d (UgriZ) photometry from the Wide Fiel d 
Survey (WFS. iMcMahon et al.ll200ll . Ilrwin fc Lewis! 1200 ll ). 
In the Lockman Hole we have 3-band photometry (gri) from 
the SWIRE photometry programme, with some additional 
U-band photometry. In Chandra Deep Field South (CDFS) 
we have 3-band (gri) photometry, also from the SWIRE 
photometry programme. Good optical data exists in other 
SWIRE fields, notably the CFHTLS^ data in XMM-LSS and 
data from ESQ W FI surveys in Elais-Sl (|Berta et al.ll2006l . 
iBerta et alj|200sh . however, these fields have not been used 
for much of our analysis as it is known that our photometric 
reds hift estimat e s are not as good in these fields. 

iFrost et al.l (|2009T ) have estimated completeness limits 
of u = 22.7, g = 23.8, r = 23.2, i = 22.4, z = 21.2 Vega 
magnitudes and /3.6um = 10, /.i.sum = 15mJy and we use 
those estimates where required in our following analysis. 



2.3 Photometric redshifts 

We use the ph o tomet ric redshifts given in 



iRowan-Robinson et al 



2008h 



These use the code ImpZ 
( Babbedge et al.l 12004 IRowan-Robinson! 120031 with up- 



dates as described in IRowan-Robinson et alj 120051 and 
iBabbedge et~ai1 l2006h which has been extensively tested 
and applied to SWIRE data. 

ImpZ is a template-fitting code that utilises the opti- 
cal and IRAC 3.6um, 4.5um detections to produce reliable 
photometric redshifts for both galaxies and AGN, account- 
ing for extinction. The code has been tested against numer- 
ous spectroscopic datasets within the SWIRE fields with 
optical data down to r < 24 and spectroscopic redshifts 
z < 3 (though with most z < 1.5). For all the samples, the 
mean systematic offset between the photometric and spec- 
troscopic redshifts was found to be negligible to the preci- 
sion of the phot ometric redshifts . For e xample, the Elais-Nl 
sample used bv lBabbedge et all (|2006 | ) has a systematic off- 
set of only (Az/(1 + z)) = +0.0037. IRowan-Robinson et all 
(|2008h provide a detailed analysis of the catastrophic failure 
rate r\ and photometric accuracy fphot = ((Az/(1 + 2 )) 2 ) 
as a function of number of photometric points, nb an d, mag- 
nitude cuts and \ 2 from the fit. From their Figure 11 we 
estimate that we would find Optot ^ 5% and r\ ^ 7% with 



1 http:/ /ssc. spitzer. caltech.edu/mips/dh/ 



2 http://www.cfht.hawaii.edu/Science/CFHTLS/ 
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( 2 < 5. and riband ^ 5 and r < 24. We use these con- 



straints for the analysis presented in Section 14.11 though 
with a much brigh ter r < 23.2 limit. However, as shown in 
iFrost et all (|20Qgh and discussed in Section 12.71 these con- 
straints are difficult to model, so we modify them slightly 
for the analysis in Section 14.21 The \ 2 constraint does not 
affect many galaxies and is difficult to adapt so we maintain 
this. We create a modified n band measure being the num- 
ber of bands that exceed the specific completeness thresh- 
olds u = 22.7, g = 23.8, r = 23.2, i = 22 A, z = 21.2 and 
/3.6u.m = 10, /4.5|im = 15mJy (as opposed to the more 
general requirement of a detection). We adopted a selec- 
tion of n(, and 4 which provided a similar set of galaxies to 
"-band ^ 5 and thus, presumably, a similar redshift accuracy. 



2.4 Optical classes and stellar mass estimation 

We also used the spectral energy distribution (SED) 
classifications and s t ellar masses given in (RR08: 
iRowan-Robinson et al.l I2008T ). They adopted a two- 
pass approach in order to fit the photometric redshift, 
optical/near-infrared SED and FIR SED making maximum 
use of the near/mid-IR Spitzer data. The optical and IRAC 
3.6 & 4.5/xm bands are fit first with a range of optical-near 
IR SEDs based on stellar population synthesis models. The 
main purpose of this first pass is to determine the level (if 
any) of excess emission in the Spitzer bands from dust, as 
well as separate AGN and galaxy spectral types for the 
second pass. The second pass includes a refitting of the 
optical SEDs to a finer redshift grid, as well as a fit to the 
far-IR component of the data with a range of mid-to-far IR 
templates. In this work we do not make use of their FIR 
luminosities or classifications but we note that because of 
this two-pass method the optical luminosity estimates and 
hence stellar masses should not be biased by any residual 
FIR contamination in the mid-IR Spitzer bands. 

The stellar masses use stellar synthesis templates (see 
Section 3, Figure 1 an d Table 2 of RR08). Start ing with em- 
pirical templates from lYoshii fe Takaharal d 19881') for galaxies 
of typ e E, Sab, Sbc, Scd, Sdm, and from lCalzetti fc Kinnevl 
(|l992T ) for starbursts, RR08 used spectroscopic data for 5976 
galaxies, for many of which they hav e 10-band photometr y 
from the CFH12K-VIRMOS survey (|Le Fevre et al.ll2004h . 
to improve these empirical templates. The latter were then 
regenerated to higher resolution using simple stellar popula- 
tions, each weighted by a different SFR and extinguished by 
a different amou nt of dust, Ay. This pro cedure, based on the 
synthesis code of lPoggianti et all |200ll ). gave the templates 
a physical validity. Minimization was based on the Adaptive 
Simulated Annealing algorithm. Detail s on this algorithm 

(|2004h . 



and on the fitting technique are given in lBerta et al 



For each galaxy they estimate the rest-frame 3.6 /im lu- 
minosity, ^L„(3.6), in units of Lq, and using the stellar syn- 
thesis models estimated the ratio (M* /Mq)/(vL v (3.6) /Lq) 
to be 38.4, 40.8, 27.6, 35.3, 18.7, 26.7, for types E, Sab, 
Sbc, Scd, Sdm, sb, respectively. [Note: measuring the 3.6 
]Um monochromatic luminosity in total solar units, not in 
units of the sun's monochromatic 3.6 fim luminosity.] Al- 
ternative estimates of M* using the B-band luminosity 
agree with our preferred method to within 10-20%. Es- 
timates based on 3.6 fim should be more reliable, since 
there is a better sampling of lower mass stars and less sus- 



ceptibility to recently formed massive stars. These mass 
estimates would be strictly valid only for low redshift. 
For higher redshifts the mass-to-light estimates will be 
lower since for the oldest stellar populations, M /L varies 
strongly with age jBruzual A. fc Charlotl 1 19931 . see their 
Figure 3). This c an be approximately modeled using the 
iBerta et al] \2QQ4 ) synthesis fits described above, with an 
accuracy of 10%, as (M*/M Q )/(vL„ (3.6) /£©)(*) = 50/[a + 



1.17(i/t ) 



where to is the present epoch and a 



0.15,0.08,0.61,0.26,1.44,0.70 for SED types E, Sab, Sbc, 
Scd, Sdm and sb, respectively. 

This approach should correctly capture the different 
evolutionary behaviour of stellar masses in star-forming 
galaxies and of early-type galaxies. 

For each galaxy we have the best fit photo-2, the tem- 
plate classification and stellar mass. There are 15 num- 
bered optical classifications Ellipticals (1-2), Sab (3), Sbc(5), 
Scd(7), Sdm (9), star-burst (11), AGN (13-15) with types 4, 
6, 8, 10 intermediate between the other optical templates. 
For much of our analysis we exclude the galaxies best fit by 
AGN templates as their photo-z are poor and the estimates 
of their star-formation and stellar mass will be strongly con- 
taminated by the AGN. 



2.5 Masks 

We have constructed a conserv ative mask b ased on those 
used in the clustering analysis of lFrost et al] (|2009l fl. Those 
clustering masks exclude regions where (1) the optical data 
is insufficient for reliable photo-z determinations (2) the 
IRAC completeness was low and/or variable and (3) fore- 
ground stars were located. Specifically we include only those 
optical frames where the rgs depth was recorded as 23.2 or 
better. This rgs depth is an estimate of the r-band 95% 
completeness limit and = 23.2 is the expected depth 
for good photometric nights. This means we are selecting 
good quality optical fields. Our optical catalogues are also 
restricted to r = 23.2. We reject areas where the IRAC 
coverage was less than 4 pointings, ensuring the complete- 
ness is above 90% and variations are less than 2% at 10/xJy. 
We also apply a coverage cut of 30 at 24 u.m giving com- 
pleteness > 90% and variation of a couple of percent for 
24 |im sources above 400 iijy. To ma sk stars we follow 
the method of IWaddington et all l|2007h in which 2MASS 
(|Skrutskie et al] 120061 )" point source catalogue sources with 
K < 12 are identified to be stars and a circular mask of 
radius, R, with log 10 (_R/") = 3.1 — 0.16K is applied. We ad- 
ditionally exclude regions where the MIPS coverage is poor 
relative to the majority of the data with a threshold chosen 
by examining the histogram of the coverage. For the MIPS 
70(J.m we use a coverage threshold of 12, while for the 160u.m 
we use a coverage threshold of 3. 

Although parts of our analysis could be carried out with 
less conservative masks we choose to apply the same mask 
to all the work in this paper so we are always comparing 
sources in the same area of sky. With our aggressive masking 
the unmasked areas of the four fields Elais-Nl, Elais-N2, 
Lockman and CDFS are Q = 4.33, 2.41, 2.75 and 1.85 sq. 
deg., totalling 11.33 sq. deg.. 



http:/ /astronomy.susx.ac.uk/~ mif20/masks/indcx.html 
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2.6 Sub-sample selection 

For part of our analysis we divide the sample into stel- 
lar mass and redshift (M* , z) cells over a range ^ z ^ 
2.0 and 9.0 < log M„/M < 12. Cell sizes were se- 
lected to provide relatively uniform number of sources in 
each cell and thus a reasonable balance between signal-to- 
noise and resolution on a scale that is easy to compare 
with other data and models. Our cell boundaries are z = 
0., 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.25, 1.5, 2 and log M»/M = 
9,9.25,9.5,9.75, 10., 10.25, 10.5, 10.75, 11,11.4, 12. The stel- 
lar mass and redshift distribution of the cells is shown in 
Figure [1] 

With redshift slices of Az « 0.2 and area of Q = 
11.33 sq. deg. our sample covers a co-moving volume 
of 6.3x 10 6 /t~ 3 Mpc 3 at z — 1. Numerical simulations 
l|Mo fc White! 120021 ) predict that this volume is sufficient 
to include the progenitors of 6 of today's 10 15 M Q clusters, 
which have a co-moving number density of 10 -6 /i 3 Mpc -3 , 
arguing that we probe a fair sample of the Universe. 



2.7 Selection effects 

Our sample and photo- z selection criteria are complex so we 
need to be cautious about selection effects. In summary, our 
selections for the stellar mass analysis are: /36(im > 10/xJy, 
r < 23.2, n' band > 4, \ 2 < 5, template classification ji < 13 
and lying outside the mask defined in Section \2. 5 1 

We are primarily investigating the mean FIR luminos- 
ity to stellar mass ratio, i.e. the specific FIR luminosity, 
and are not concerned with the total number of sources or 
the luminosity density. Thus in any M» , z cell we only need 
to worry about selection effects that affect sub-populations 
within the cell in different ways. In other words our estimate 
of the mean specific FIR luminosity in a cell would be unaf- 
fected if we randomly exclude 50% of all galaxies in the cell 
but would be affected if we exclude 50% of the most FIR 
luminous and none of the others. 

The x 2 selection might affect some galaxies types more 
than others, however, with a broad range of templates, care- 
fully modified to fit the SWIRE populations, we expect this 
cut will have a smaller effect than the flux related criteria. 

The number of bands criteria and the flux cuts do need 
to be considered. Since different galaxies have different K- 
co rrections this introd uces some differential effects. As noted 
bv lFrost et al.1 ll2009h . in a cell of fixed stellar mass (or 3.6 
H-m luminosity) and redshift the optical cuts will exclude 
redder galaxies. As these galaxies typically have lower spe- 
cific star-formation rates ignoring this effect would artifi- 
cially inflate our specific FIR luminosity estimates. To ac- 
count for this we weight the galaxies to account for this 
incompleteness and flag as incomplete cells where our in- 
completeness corrections may be inadequate. 

For all cells we use a V max weighting. The model tem- 
plates give us a K-correction in every band. For each galaxy 
we use our measured K-correction and observed magnitude 
to calculate the maximum observable redshift z max ,i in each 
band, i. We then compute z max taking into account our 
3.6um, r-band and n' h!lnd constraints. The available volume 
can then be calculated taking into account the cell limits 
and z max . We also flag as incomplete any cells in which for 
any template; (1) we found no examples of galaxies best 



fit by that template and (2) a galaxy with that template 
could not, in principle, have been found in the cell, given 
the quoted selection limits. Later, we also exclude any cells 
with fewer than 100 galaxies. 



3 STACKING METHOD 

Stacking analysis co-adds the signal in a map at the po- 
sition of a class of galaxies, allowing for the reduction of 
noise associated with measurements of individual g alaxies, 
whether this be confusion noise (e.g. ICondonlll974l ) or in- 
strumental/background noise. The level of noise reduction is 
governed by a number of factors, including the ability to cat- 
egorise the target catalogues into groups with similar proper- 
ties and the resulting number of targets on which the stack is 
perf ormed. Stacking has been used at sub-mm wavelengths 
(e.g.lPeacock et alj200d : iDve et al.ll2007l ; iTakagi et af]|2007l ; 
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Spitzer data ( Dole et al.ll2006l ). Stacking has also been used 
to investigate the far-infrared spectral energy distributions 
of Spi tzer galaxies in sm all fields at z ~ 0.7 (jZheng et al.l 



l2007bl l and 1.5 < z < 2.5 (|Papovich et alj|2007 

The conventional stacking technique is to extract a re- 
gion from the map around each target source and stack these 
together to produce an average image. A mean background 
is subtracted; usually estimated either from the global map 
or from the extremities of the average image. The average 
flux is then determined by either calculating the total flux in 
some aperture or fitting the point-spread function to the im- 
age. If the noise in the original maps is uncorrelated with the 
target sources (being either instrumental noise, background 
thermal noise, confusion noise from unrelated galaxies, or 
foregrounds such as Galactic cirrus) then the noise in the 
average map will be reduced by a factor of y/W where N 
is the number of targets. Furthermore, if N is large then 
the noise will approach a Gaussian distribution due to the 
central limit theorem. Under the assumption that the tar- 
gets are isolated and point-like at the telescope resolution 
a suitably weighted point-source profile fit is the optimal 
estimator of the average flux of the targets. 

An additional contribution to the average flux will come 
from sources that are spatially correlated with the target 
galaxies. This is a difficult contribution to model as it de- 
pends on the correlation function of galaxies which may be 
luminosity or type dependent. We discuss this further in 
Section f3. II In order to minimize this bias we fit the point- 
source profile and a constant background simultaneously 
over a limited radius. (A similar technique for source detec- 
tion and photometry has been discussed bv lSavage fc Oliver! 
l2007f ). The simultaneous fit means that the background is 
estimated from a local area that includes the region under 
the source, in contrast to methods using a sky annulus which 
exclude the source region. The point-source fitting and lim- 
ited radius means that the central parts of the image, where 
the source flux dominates any correlated background, have 
a greater weight. 

We performed the fit in ID using a radial point-source 
profile P(r) estimated from the 2D calibration files. The two 
parameter intensity profile function I(r) — Io P{r)/ Po + b 
was fit, limiting the data to that within the first Airy disk 
minimum. We used a minimum x 2 fit with errors in the 
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Figure 1. Stellar mass vs photo-z plane for groups in our analysis. Each point represents the centre of a M t ,z cell with the size 
representing the number of objects in the cell according to the code in the legend. 



mean image intensity estimated from the scatter between 
the individual images in the stack. This would be optimal 
if the population variation is small and the errors are un- 
correlated, so it probably under- weights the central pixels 
but is adequate for our purposes. From Io we deduce the flux, 
/, using the effective beam calibration factor, fibeam, given 
earlier. Looking at the x 2 we find reasonable fits whenever 
we have 9 or more galaxies. 

It is simple to use the central limiting behaviour of the 
stacking technique to get a good estimate of statistical er- 
ror in the flux. However, we are concerned about systematic 
noise terms and so will use an error calculation based on 
the variation from field to field (see Section [A. Our tech- 
nique can easily be extended to include correlated errors, 
to accommodate a model for any background (e.g. one from 
correlated sources) and also to provide an estimate for the 
variation in the population flux as well as the mean. 

3.1 Simulations 

From the beam sizes calculated above we see that the density 
of beams on the sky is 21000 and 4700 per square degree 
at 70um and 160 um, respectively. Comparing these with 
Table [1] we can see the number of sources per beam can be 
high, particularly at 160u.m. We should thus be concerned 
about confusion i.e. where correlated neighbouring sources 
spuriously increase the stacked flux. 

Our method aims to mitigate this problem by using 



the simultaneous source and background fitting. To test this 
we have run some simulations. We need to investigate the 
behaviour in high and low source densities regimes and to 
include galaxy clustering. Our approach is to use the real 
catalogue positions. For the low source density case we take 
the sources with S24 > 400/iJy. We are not interested in the 
absolute fluxes, so we model the long-wavelength fluxes as 
SVo = Sieo = <S , 24- We insert these fluxes at the catalogue 
positions and convolve with the corresponding 2D point re- 
sponse function. For the high source density case we repeat 
this using sources with S3.6 > 10/xJy. We then undertake 
our stacking analysis using the same stellar mass and red- 
shift cells. We do not add noise to the simulations as we are 
concerned here with systematic biases from confusion. The 
results are shown in Figure [2] 

These simulations show that the uncertainties increase 
as the mean flux in the stack decreases. The results are 
roughly similar for the two different sample densities and 
the two different beam sizes. At any given flux the uncer- 
tainties appear to be larger for smaller stack samples. If we 
only consider the larger stack samples then there is some in- 
dication of systematic underestimation of the simulated flux 
at faint fluxes. This effect appears to be less than 0.5 dex. 
At first glance this is surprising since our naive expectation 
was that the correlated signal would increase the fluxes. This 
suggests that our background estimation is biased upwards 
by the correlated flux. In principle we could model this bias. 
However, our simulation is limited as (1) we have assumed 
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a direct correlation between the long wavelength fluxes and 
short wavelength fluxes and (2) we have ignored the con- 
tribution from sources with higher space densities than the 
SWIRE 3.6 urn sample. The first assumption is difficult to 
model as it relies on an understanding of both the luminosity 
and clustering properties of the FIR emitting galaxies which 
are poorly constrained. As we show in Table [T] the second 
assumption is relative modest as we appear to resolve much 
of the FIR background within the SWIRE sample. However, 
we note that the SWIRE catalogues will not have resolved 
sources closer than the 3.6um beam and so will underesti- 
mate the number of close pairs. Had such close pairs been 
included it would have counteracted the observed bias. 

We conclude from the simulations that systematic ef- 
fects due to source correlation in highly confused regions 
may underestimate the fluxes at low fluxes, but by a factor 
smaller than 0.5 dex.. 



3.2 From stack fluxes to background intensities 
and luminosities 

The stacking technique naturally gives us the mean flux of 
galaxies in a class. 

To derive the contribution of that class to the back- 
ground intensity (ul v ) we take the flux per unit frequency 
and divide by the unmasked area of the catalogue from 
which that class was selected i.e. vl v = vy~] f^/Q. 

The conversion to specific luminosity is more in- 
volved. To minimise differential effects across the cell and 
to include the appropriate completeness corrections we 
multiply each image in the stack (arising from galaxy 
i) by a scaling factor, u>i, before stacking (with Wi = 
AnDLizifKFmiz^V^^M^). Then, rather than dividing 
the total stack by the number of galaxies, we divide by 
X^Knaxi- The lu mmos ity distance, Dl, and K-correction, 
^fir, transform flux into FIR luminosity while the stel- 
lar mass, M*i, converts to specific luminosity and the Vmax 
terms correct for the incomplete sampling of the cell volume. 

To calculate -Kfir we use a mean spectral energy dis- 
tribution (SED) averaged over all galaxies in the cell. In 
Figure [3] we compare t he 70/160 colo u rs of all the cells with 
model templates from iPolletta et al.l (|2007l ) and individual 
galaxies detected in both bands. The individual detections 
have similar colour distribution to the stack samples but 
with a tail to warmer colours. This tail is expected because 
the detection criteria selects the most FIR luminous objects 
and there is a well known correl ation between FIR lum inos- 
ity and dust temperature (e.g. IChapman et alJl2003l ). No 
single template fits all the samples at all redshifts. How- 
ever, we see that the Sc template (plotted with a thicker 
line style) provides as good a fit as any other template over 
the range of classes and redshifts. This is natural as we'd 
expect the mean SED to be roughly the same as a galaxy 
with a moderate level of FIR activity. 

We thus adopt an Sc template as a reasonable compro- 
mise between the different SEDs and use this to compute 
the FIR luminosity, Lpm, where we define the FIR range 
to be from 5um ^ A ^ lOOOum, from 70 or 160 urn fluxes. 
It should be noted that the Sc template peaks at lOOum 
(in vf v ) or 130u.m (in fu) so the K-corrections at the two 
different wavelengths are different (in opposite directions 
for z < 0.2). Figure [4] shows a wide variety of empirically 



based spectral energy distributions for galaxies detected in 
infrared bands, we see that the 160um band on its own is 
reasonable bolometric power indicator, i.e. any uncertainties 
in which template is appropriate have a small effect on the 
K-correction. 

If the far infrared luminosity traces the star-formation 
activity then the specific FIR luminosity is a measure of the 
specific star formation rate. However, the FIR luminosity 
also traces emission from AGN and from diffuse dust heated 
by the ambient stellar radiation field ("cirrus"). The AGN 
emission is warmer than either star-formation or cirrus and 
is expected to be less significant at these long wavelengths. 
This is confirmed by our SED analysis in Section 14.11 The 
well-known FIR- Radio correlation (e.g. ICondonlll992f ) sug- 
gests that even the cool "cirrus" emission indirectly traces 
star-formation. Strong AGN with modest extinction are ex- 
cluded by the optical SED modeling. Weaker or more ob- 
scured AGN will remain and provide some level of contam- 
ination. The level of this contamination is hard to estimate 
but typically AGN fractions are found to be ~ 30% (e.g. 
IPolletta et al.l liood ) . A strong AGN contamination would 
reveal itself in differences between our estimates arising from 
70 and 160 urn. Thus it is plausible to relate the total FIR 
luminosity to the obscured star formation rate. We relate 
the FIR luminosity to the total star formation rates (i.e. 
including an estimate o f the unobscured contribution ) us- 
ing the conversion from Rowan-Robinson et al.l (1 19971 ) and 
iRowan-Robinson et al.l (|2008l ) expressed as: 



^ = 0.51 x 10 10 



SFR 



AL 



©yr" 



where we have taken the fraction of UV energy absorbed 
by dust to be e = 2/3, and Z/fir/£60 = 1.67. These values 
are appropriate for an M82 like spectrum. It is worth not- 
ing that if the obscuration is less than this, as you might 
expect for systems with less active star formation, then we 
will underestimate the star formation rate. 



3.3 Error estimation 

To estimate the systematic uncertainties we calculate our 
errors using the variation in mean fluxes we get from field 
to field. This approach accounts for sampling variance er- 
rors and some of the systematic errors arising from our use 
of photometric redshifts. Although our photometric redshift 
technique is the same across all the fields, the optical data 
comes from different telescopes; this means the field to field 
variations could be significant if the photo- z have any subtle 
dependencies on the bands or optical limits. 

The areas of each field fij are given in Section 12.51 
For each sub-sample we calculate an average flux (or 
specific star-formation) weighted by these areas / = 
y\ Clj fi/y)- We estimate the error on the resulting 
weighted average as 



The first term is a factor to scale from the population 
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Figure 2. Ratio of the output to input flux ratios for simulations of our stacked samples. Flux units are arbitrary. We use the stack 
samples from the Elais-Nl field (our largest field, the total sample has over twice this area). The colours of the symbols indicate whether 
we placed the FIR fluxes at the location of 3.6 (im or 24 um sources and whether we used the 70 or 160 urn beam. Symbol sizes relate 
to sample sizes in the same way as Figure [T1 except the symbol sizes are reduced by a factor of two. Cells with fewer than 10 sources are 
marked with crosses. 



variance to the variance in the mean, the next two terms are 
the weighted estimator for the population varianceQ 

Any stacks with fewer than 9 galaxies (the threshold for 
reasonable \ 2 found in Section [3| were excluded from the 
average and any cells that had fewer than 100 galaxies after 
the averaging were excluded. 



4 RESULTS 

4.1 Contributions to CIRB from various 
observational sub-samples 

Before we examine the specific far infrared luminosity we 
explore the contribution of different subsets of the SWIRE 
catalogues to the cosmic infrared background (CIRB). These 
are tabulated in Table [T] 

The SWIRE galaxies which are bright enough to be 
detected individually at 70 or 160um (rows 3 & 4) contain 
only a small fraction (< 12%) of the CIRB information of the 
SWIRE data set as a whole (row 1) which provides strong 
motivation for our stacking analysis. The optical sub-sample 



4 e.g. http://pygsl.sourccforgc.net/rcfcrcncc/pygsl/nodc36.html 



(row 5) contains about 60% of the FIR information. This se- 
lection is similar to that used in most of this paper. It is en- 
couraging that we detect this much of the FIR flux, however, 
it emphasises the need for deeper optical data to fully ex- 
ploit the SWIRE data. The SWIRE 24 um catalogues detect 
about half of the CIRB seen in the 3.6uxn sources. We find 
that S WIRE 24um catalo gues resolve 30-40% of the CIRB 
seen by iDole et al.l (|2006t ) in fainter 24um samples, while 
the SWIRE 3.6um catalogues resolve 70-80%. Estimation of 
the total CIRB from direct measuremen t is highly uncer- 
tain but if we use the estimates quoted bv lDole et al.l (|2006t ) 
then the SWIRE 3.6um catalogues resolve about 70% of 
the 70um background and 60% of the 160um background. 
Roughly half of the background is resolved in the optical 
samples at r < 23.5. 

Some other points to note from this table are that the 
FIR colours generally become warmer as we move from early 
to late types, starbursts and AGN (as expected). As we 
noted before when discussing Figure[3l the colours of sources 
catalogued at 70 or 160uxn are warmer than those selected 
at other bands. This is because these sources tend to have 
higher FIR luminosity and are thus preferentially star-burst 
galaxies or AGN. This effect is less pronounced at 160um 
which picks up cooler sources. 
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Figure 3. f v colours ratios at 70/160 uin for stack sub-samples. Filled circle data points arc stack samples, colour-coded according to 
stellar mass using the same scheme as for Figu res [9l 1111 Black data points are catalogued galaxies extracted from the SWIRE Elais Nl 
field. Model templates are taken from lPolletta et al. I 1120071 1. 



Selection 


Number 


Number Density 


ho 


1 160 


ho/heo 






/ 1000(sq. deg.)- 1 


/nWm _2 sr _1 


/nWm~ 2 sr -1 




536 > 10/iJy 


394014 


34.8 ±1.6 


4.231 ±0.085 


8.77 ±0.44 


0.48 ±0.03 


S 24 > 400AtJy 


21146 


1.87 ±0.05 


2.159 ±0.066 


3.17 ±0.06 


0.68 ±0.02 


S 7 o > 30mJy 


856 


0.076 ±0.004 


0.703 ±0.039 


0.68 ±0.03 


1.04 ±0.08 


5i60 > 90mJy 


877 


0.077 ±0.003 


0.542 ±0.041 


0.71 ±0.04 


0.76 ±0.07 


S 3 e > 10/iJy; 13.5 < r < 23.5 


223402 


19.7 ±0.6 


3.247 ±0.099 


5.71 ±0.32 


0.57 ±0.04 


13.5 ^r< 23.5; X 2 < 5; riband > 4 


182627 


16.1 ±1.0 


2.390 ±0.235 


4.28 ±0.38 


0.56 ±0.07 


E (j 2 = 1,2) 


29700 


2.62 ±0.09 


0.260 ±0.023 


0.56 ±0.05 


0.47 ±0.06 


Sab (j 2 = 3, 4) 


9727 


0.86 ±0.05 


0.198 ±0.013 


0.38 ±0.02 


0.52 ±0.05 


Sbc (j 2 = 5, 6) 


23306 


2.06 ±0.12 


0.417 ±0.043 


0.80 ±0.06 


0.52 ±0.07 


Scd (j 2 = 7,8) 


54864 


4.84 ±0.20 


0.911 ±0.070 


1.63 ±0.12 


0.56 ±0.06 


Sdm (j2 = 9, 10) 


37594 


3.32 ±0.28 


0.293 ±0.056 


0.45 ±0.08 


0.65 ±0.17 


Star-burst (j 2 = 11) 


23549 


2.08 ±0.17 


0.276 ±0.078 


0.46 ±0.08 


0.60 ±0.20 


AGN {j 2 = 13 - 15) 


3887 


0.34 ±0.12 


0.053 ±0.019 


0.06 ±0.03 


0.89 ±0.51 


S 24 > 60/iJy 






5.94 ±1.02 


10.72 ± 2.28 




CIRB 






6.4 


15.5 





Table 1. Contributions to the 70nm and 160um FIR background from various populations. All catalogues and maps have been masked 
by the mask described in the text. The first four rows are galaxies extracted from the Spitzer data only catalogues. The fifth row comes 
from the Spitzer /optic a l cros s-matched catalogues. The remaining rows above the line are extracted from the photo-z catalogues of 
iRowan -Robinso n et al. | j- 1 ' ">-| > which has been filtered to remove galaxies with poor quality photo-2 and classifications (i.e. those with 
n band < 4, or \ 2 ^ 5). The data are averaged over four fields Elais-Nl, Elais-N2, Lockman, CDFS weighted by the unmasked areas of 
4.33, 2.41, 2.75, 1.85 sq. deg. respectively for a total of 11.33 s q. deg. Errors are deduced from the field to field variations. The last two 
rows are estimates extracted from Table 1 of lDole et al. (2006) 
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Figure 4. Bolometrically normalised spectral energy distributions (SEDS) for a range of models from lXu etal] Jl998h . SEDs are shown 
in the rest-frame and are normalised to have the same power over the range 5 < A/u.m < 1000. The models include starburst galaxies, 
normal spirals and AGN. Overplotted are the 3 Spitzer MIPS bands (24, 70 & 160|J.m) for a galaxy observed at z = 1. Notice that the 
bolometric normalisation is similar to normalisation at 100|J.m. Our analysis uses 70 and 160|J.m, much previous work uses 24|j.m. 



4.2 Specific far-infrared luminosity as a function 
of stellar mass 

To explore the specific star-formation we measure the ratio 
of far-infrared luminosity to stellar-mass i.e. the specific FIR 
luminosity. An early exploration of this as a function of op- 
tical luminosity ( as a proxy for stellar m ass) from ISO data 
was presented bv lOliver fc Pozzll (|2005f ). We use the stellar 
mass and redshift cells described in Section \2. 61 

In Figures [5] and [6] we plot the specific FIR luminosity 
as a function of stellar mass. We omit points where the frac- 
tional error (calculated using the field-to-field variations) is 
more than 1. Data points flagged as incomplete (Section l2,7[) 
are unfilled. 

We see similar results at both 70 and 160 \im. This 
agreement is an important validation of our method as the 
data are completely independent and the K-corrections ap- 
plied are very different. We are thus confident in combining 
the two data sets to give a single estimate of the specific 
FIR luminosity (being a weighted average of the estimate at 
each wavelength). The scatter between the two independent 
wavelengths is indicated in the error bars. 

In Figure [5] we compare our average specific FIR- 
luminosity in our lowest redshift bin (where we cover the 
widest mass range) with other data from the literature. We 
see a simple power-law trend but a shallower slope, particu- 



larl y with higher spec ific star-formations at higher masses. 
The lChen et all (|2009l ) work models the SDSS spectra and is 
sensitive to star-formation over a longer time-scale (~ lGyr) 
than the FIR. This may lead them to find higher specific star 
formation rates due to the evolution over that time-scale. 
Alternatively, this may be an obscuration effect, if the op- 
tical observations miss heavily obscured star-formation and 
if that deficit is greater in higher mass systems. This might 
be expected if star formation in massive systems tends to be 
in more d eeply embedd e d site s as is seen in Arp 220. The 
analysis of lDamen et all i|2009l ) includes star formation mea- 
sures from 24 urn and so should include obscured star for- 
mation, but may miss cooler contributions that we pick up 
at longer wavelengths. This latter explanation i s tentatively 
suppor ted by the IRAS measurements from IZheng et all 
(|2007bf ) at lo nger wavelengths whi ch lies between our results 
and those of lDamen et al. | (I2009T). The remainin g discrep- 
ancy between our work and Zheng et all (|2007bl ) (at lower 
z) could then be explained as due to evolution across our 
redshift bin, which will preferentially increase the sSFR in 
the higher mass bins. A final possibility is that there are 
inconsistencies in the stellar mass estimates, either through 
modelling or through photo- z estimates. This is likely to be 
more of a problem for us at these lower z, where photo-z 
errors are more significant, than it is at higher z. 

Comparing with data at higher redshifts (Figure (6]l the 
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Figure 5. Specific FIR luminosity or specific star-formation rate as a function of stellar mass in redshift classes. Rcdshift ranges from 
< z < 0.2. Filled circles are the average of our estimates from the 70 fcl60|xm data which are the weighted average over four fields 
(Elais Nl, Elais N2, Lockman and CDFS). Thick vertical error bars come from the field-to-field scatter and thin error bars with hats 
from the variation between 70 & 160|J.m estimates. All points have used at least 100 galaxies in the stacking analysis. Solid lines are our 
power-law trends with parameters give in Table [2] Calibration to specific star-for mation (right-had axi s) is given in the text. Previous 
estimat es of specific star- fo rmations rates are estimated from 24 |xm data (star s) iDamen et alj (1200911 (0.1 < z < 0.3), IRAS (open 
circles) IZheng et al. | d2007aT) z < 0.02, and from SPSS spectra (triangles) z ~ O.l lChen et alj ll2009l )7 Dashed lines are predictions from 
the semi-analytic models of lDe Lucia fc Bl aizot (200]]). 



agreements are much better. Our large area means we are 
able to divide our data into finer redshift bins than previous 
work while still maintaining high statistical precision. We 
find similar slope and a mplitude to obscured tracers e.g. esti- 
mated from 24 urn data (|Zheng et al.ll2007bh . Their observa- 
tions were over a much smaller field and so they were limited 
to lower ste llar masses. A t z > 1 we find sSFR significantly 
higher than lDamen et all (|2009l ) ( who included FIR i ndica- 
tors) , but similar to those found bv lDunne et al] (l2009t) (who 
use ra dio). Our slope is hard to compare with Dunne et all 
|2009l) as their redshift ranges are much larger and differen- 
tial ev olution across their bin will be significant. [Elb az et alj 
(|2007l ) derived an empirical relation from 24um and UV 
data and quote SFR/(M Q yr _1 ) = 7.2[Af,/(10 10 M Q )] ' 9 at 
0.8 < z < 1.2, ag ain limited to lower stellar masses. The 
[Elb az et al.1 (|2007r ) measurement has a consistent amplitude 
in the limited range where our stellar mass ranges overlap. 
However, th e y have a shallower slope than ours and that of 
IZheng et al] (|2007bf) ; this might be an effect of differential 
evolution across their wide redshift bin. 



for consistent measurements over a wide dynamic range in 
mass and in narrow redshift bins. 

There is a strong evolutionary trend which we explore in 
more detail in Section 14.31 Looking within each redshift bin 
(while bearing in mind that we only see lower mass objects 
at lower redshifts) we see that there is a consistent decline 
in specific FIR luminosity as we move to higher stellar mass. 
We model this trend as 



sSFR^^oc 



M 



1O X1 M 



This power-law model is a reasonable fit and parameters for 
the average fit are given in Table [2] There is some variation 
in the slope with redshift, the slope is steeper for 0.6 < 
z < 0.8, but a mean value of /3 = —0.38 ± 0.14 provides a 
plausible description of the data. 

For a direct comp arison with semi-ana l ytic m odels we 
take the simulations of |Pe Lucia fc Blaizotl (|2007l fl and se- 
lect galaxies in the same stellar mass and redshift cells as 



Overall it seems that the relations are smooth and there 
is no strong evidence for any breaks. However, there is a need 



5 Publi cly available on the Millennium Simu lation data download 
site (see lLemson fc Virgo Consortiumll2006l) . 
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Figure 6. Specific FIR luminosity or specific star-formation rate as a function of stellar mass in redshift classes. Mean points are the 
weighted average over four fields (Elais Nl. Elais N2, Lockman and CDFS) and over 70 &160u.m data. Error bars are estimated from 
the scatter between fields and bands as in Figure [5] Incomplete cells where not all galaxy templates are represented are indicated with 
open circles. Power law fits to complete data are shown with solid lines. Comparison samples from similar redshift ranges are plotted 
with fainter symbols. Since the redshift bins are not identical these sometime appear on more than one plot or with data from two 
bins in the original work on the same plot. E stimates for combined 24um and GALEX UV photo metry come fromlElbaz et al. I J2007t l 
(shaded, 0.8 < z < 1.3) IZheng et all j2007al ') (circles, z bins, 0.2-0.4, 0.4-0.6, 0.6-0.8,0.8-1.0) and iDamen et al.l j2009h (stars, z bins 
0.1-0.3, 0.3-0.5, 0.5-0.7, 0.7-0.9, 0.9-1.1, 1.1-1.3, 1.3-1.5, 1.5-1. 7, 1.7-1.9 ) Radio estimates from lDunne et al.l j2009h (down triangles, z 
bins 0.2-0.7, 0.7-1.2, 1.2-1.7, 1.7-2.2) and from DEEP 2 spectra lchen et al.l i2009l ) (up triangles, z ~ 1). The parameters of the pow er-law 
fits to are data are tabulated in Table [2] Dashed lines are predictions from the semi-analytic models of |Pe Lucia fc Blaizotl J2007t) 



our data. This model includes AGN feedback which effec- 
tively suppresses star formation for galaxies with large cen- 
tral black holes (typically galaxies with M> 10 10 ' 5 M©). The 
comparison is shown in Figure [S] for our lowest z bin and 
in Figure [6] for the others. This model fits our low mass 
data well, but the model under-predicts our high-mass end. 
The same model agrees better with other data sets at the 
high mass end but not as well at lower mass. The shape of 
the model curves are identical at higher redshifts with just 
the normalisation increasing. The clear break in the model 
curves arises at the mass scale where the AGN feedback be- 
comes important. We don't see such a break in any redshift 
bin. In addition, this model fails to predict the amplitude of 
the evolution with redshift which we discuss in Section T4. 31 

We perform the same analysis for galaxies separated by 
the optical/NIR SED class. We show the separation into 
late-type galaxies (SED types 3-11) in Figure [7| and early- 
type galaxies (SED types 1-2, Figure [S| with fit parame- 
ters in Table [2] It is immediately striking that the relation 
ships are much steeper (though naturally lower amplitude) 



for the early type galaxies /? ~ —0.46 q.v. /? ~ —0.15 for 
the late-type galaxies. This may indicate that early type 
galaxies harbour some components of star-formation that 
are unrelated to the dominant host galaxy (e.g. through ac- 
cretion of gas rich companions). However, this cannot be the 
complete picture as totally uncorrelated star-formation onto 
dead hosts would produce a steeper slope with p = — 1. 

We have performed the same analysis on all of the dif- 
ferent optical galaxy classifications and these show a logical 
progression between the behaviours seen in the early and 
la te groupings shown here . A figure illustrating this is shown 
in lRoseboom et alj l|2009l ). 

If we ignore the highest redshift bin, for which the slope 
is poorly constrained, there is some indication from these 
figures and Table [2] that the slope of the relation between 
sSFR and stellar mass becomes steeper with increasing red- 
shift for both early and late type galaxies and for all galaxies 
combined. 
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Figure 7. Specific FIR luminosity or specific star-formation rate as a function of stellar mass in redshift classes for "blue" galaxies, 
i.e. those with SED types 3-11. Points are the weighted average over four fields (Elais Nl, Elais N2, Lockman and CDFS) and over 70 
&160|rm data with error bars are as for Figures [5] and [6] Incomplete points arc unfilled. The parameters of the power-law fits to are data 
are tabulated in Table [2] 
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Table 2. Power-law fits to the specific star-formation as a function of stellar mass as plotted in Figures 1 5181 Modelled as sSFR = 
Y (M /lO 11 Mq)' 3 in redshift ranges indicated. Combined 70 and 160 urn data. Averages and standard deviations of /3 are given below 
the line. Fits are calculated for all galaxies and separately for "blue" (templates Sab-Sdm and starburst, 3 ^ j% ^ 11) and "red" galaxies 
(both E templates, J2 ^ 2). Reduced \ 2 are quoted if the number of points used in the fit is more than 2. 



4.3 Specific FIR luminosity as a function of 
redshift 

By plotting the same data as a function of redshift but in 
stellar mass classes the evolutionary trend is readily appar- 
ent, see Figure [9] It is immediately clear that there is a 



dramatic increase in FIR to optical ratio (or specific FIR lu- 
minosity) by a factor of > 100 over the interval ^ z ^ 2.4. 
This is seen at all masses apart from the lowest mass bin 
which has poor statistics and a limited redshift baseline. 
Again, we see the same trends in both the 70 and 160u.ni 
data. This is an important corroboration as the FIR data 
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Figure 8. Specific FIR luminosity or specific star-formation rate as a function of stellar mass in rcdshift classes for "red" galaxies, i.e. 
those with SED types 1-2. Points are the weighted average over four fields (Elais Nl, Elais N2, Lockman and CDFS) and over 70 &160u.m 
data with error bars are as for Figures \E\ and [6] Incomplete points arc unfilled. The parameters of the power-law fits to are data are 
tabulated in Table 



are independent and the K-corrections will be different. We 
have fit this trend with a simple model 

sSFR = X(l + z) a . 

The best fits are shown in Figures O [TO] and [TT] and the 
parameters and goodness of fit are given in Table [3] 

The simple power-law fit provides a good description of 
the data in most cases. As already seen in Section [4.21 the 
zero redshift specific FIR luminosity, i.e. the parameter X, 
declines with increasing stellar mass. We see no strong indi- 
cation that q increases with increasing mass (in the range 
1O 1O ' 5 M < M < 10 12 ), as would be expected in the "down- 
sizing" scenario. 

We als o comp are with the observations from 
I Zheng et al.l (|2007bl ). Their data shows similar evolu- 
tionary behaviour at lower stellar masses. If we assume 
that the variations of evolutionary rate with stellar mass 
are statistical variations then we can make a simpler model 
with the same a for all classes. We take the average of 
the a estimates in Table rj for M > 1O 1O ' 5 M . From this 
we deduce that the Lfir/M* ratio varies as (1 + z) 4,4±0 ' 3 , 
with the error bar being the standard deviation of the mea- 
surements. The lowest mass bins are, however, inconsistent 
both with this mean evolut i onary rate and the lower mass 
estimates from IZheng et all (|2007bl V 

We show the redshift trend for the SAMs in Figure [9j 



The SAMs show specific star formation rates decreasing 
smoothly with cosmic time as the gas supplies are declining 
and the stellar masses are building up. The trend with red- 
shift in the SAMs has a mean slope of a = 2.4. which is con- 
siderably shallower than the observed data (a = 4.4 ± 0.3). 

We do the same analysis for late/blue (Figure [10|) and 
early/red galaxies (FigureQT]) as before. It is striking that for 
the late-type galaxies the variation with stellar mass seen in 
Figure [9] is considerably reduced and that there is negligible 
difference between galaxies with M > 10 10 ' 5 Mq. Assuming, 
therefore, that they are measuring the same relation, we av- 
erage the fits in Table|3]to estimate log 10 (sSSFR)/Gyr _1 = 
-1.36(±0.20) + 3.36(±0.16) log 10 (1 + z), where the error 
bars are errors in the mean calculated from the standard 
deviations given in Table [3] 



5 DISCUSSION 

Our stacking analysis shows that the specific FIR luminosity 
declines as a function of increasing stellar mass and evolves 
strongly with redshift. Before we consider the implications 
of this it is worth drawing together some of the caveats. Hav- 
ing corrected for the optical incompleteness the remaining 
selection effects are expected to be weak and should only ef- 
fect the lowest luminosity or highest redshift cells. The bias 
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Figure 9. Specific FIR luminosity or specific star-formation rate as a function of redshift in stellar mass classes. Filled circles are the 
average of our estimates from the 70 &160u.m data which are the weighted average over four fields (Elais Nl. Elais N2, Lockman and 
CDFS). Thick vertical error bars come from the field-to- field scatter and thin error bars with hats from the variation between 70 & 
160um es timates. Solid lines are ou r power-law trends with parameters give in Table [3] The shaded region represent the semi-analytic 
models of lDe Lucia fc Blaiz ot (2007) with the upper and lower bounds being the predictions for M* = 1O 1O M0 and 10 11 Mq. The dashed 
line shows the inverse Hubble time, galaxies above this line are 'bursting' producing stars at a higher rate than their historical average. 
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Table 3. Fits to specific star-formation rates as function of z in stellar mass bins. Fit is to the function sSFR = X(l + z) a . Averages 
and standard deviations for a from data with M > 10 10 5 Mq are shown below the line. Calculated for all galaxies and separately for 
"blue" (templates Sab-Sdm and starburst, 3 ^ j% ^ 11) and "red" galaxies (both E templates, J2 ^ 2). Reduced \ 2 are quoted if the 
number of points used in the fit is more than 2. a for SAM model varies between 2.44 and 2.67 over same mass range. 



against very obscured objects would act to reduce either of 
the basic trends we found. The possible bias to underesti- 
mate faint fluxes exposed in the simulations may bias us 
to slightly shallower slopes of sSSFR vs stellar mass. Sam- 
pling variance is minimised by our large survey volume and 
is included in our error bars which take into account the 
variation from field to field. Excluding galaxies with signs of 
AGN in the optical SED may bias us against star-formation 



correlated with AGN activity, or to higher AGN contami- 
nation in more massive systems where the signature of the 
AGN in the optical SED may be masked. The concordance 
between our independent analysis at 70 and 160 urn also 
argues against any of these effects being significant. Our 
biggest remaining caveat is that we are strongly reliant on 
the accuracy of the photometric redshifts and resulting stel- 
lar masses from iRowan- Robinson et~aH (|2008l ) but we have 
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Figure 10. Specific FIR luminosity or specific star-formation rate as a function of redshift in stellar mass classes for "blue" galaxies, i.e. 
those with SED types 3-11. Filled circles are the average of our estimates from the 70 &160um data which are the weighted average over 
four fields (Elais Nl, Elais N2, Lockman and CDFS). Thick vertical error bars come from the field- to-field scatter and thin error bars 
with hats from the variation between 70 & 160um estimates. Solid lines are our power-law trends with parameters give in Table [3] The 
dashed line shows the inverse Hubble time, galaxies above this line are 'bursting' producing stars at a higher rate than their historical 
average. 



mitigated against systematic errors in these by applying con- 
servative constraints on the photometric redshifts and from 
our field-to-field comparison. 

The specific FIR luminosity we measure includes ener- 
getic contributions from star formation, AGN and the am- 
bient stellar radiation field which we relate to specific star 
formation rate. Assuming these different contributions track 
each other with redshift then evolutionary trends will be un- 
affected. 

With these caveats we consider that the trends we 
have observed in specific FIR luminosity mirror trends in 
specific star formation rates. The specific star formation 
rate is dependent on the factors triggering star formation, 
the availability of gas supplies and the importance of feed- 
back processes that regulate the star formation. The de- 
cline of specific star formation rate with stellar mass (which 
has already been observed at other wavelengths) may be 
due to declining resources, i.e. more massive galaxies have 
locked more of their baryons into stars in the past. Feed- 
back also plays a role, e.g. AGN feedback in t he more mas- 
sive galaxies may suppre ss star formation (e.g. iBower et al.l 
l2006l ; ICroton et alj|2006h . However, it seems that the abrupt 
changes in specific star-formation resulting from some mod - 
els of AGN feedback (i.e. those of lDe Lucia fc BlaizotJkoOTl ) 



are not supported. So either the impact of AGN feedback 
on specific star-formation is marginal or the AGN feedback 
is not as differentially dependent on stellar mass as in these 
models. 

The evolutionary trend is both remarkably consis- 
tent across different stellar masses in the range 10.5 < 
log 10 M*/Mq < 12, and stronger than the models. This dis- 
crepancy between the models and data could arise from an 
evolution in the feedback prescription which would need to 
be less restrictive at higher redshift. Alternatively the trig- 
gering mechanisms (e.g. environmental effects) or decline in 
gas supplies in the models may need adjusting. It is worth 
noticing that enhanced evolution of the star formation rates 
enhances the change in specific star formation both through 
the star-formation directly but also through the more rapid 
stellar build-up. 

It should also be noted that the parameters in the 
models including the ones controlling the AGN feedback 
strength, were only adjusted to reproduce galaxy proper- 
ties at redshift zero and any resemblance at all to observed 
data at high redshift might be viewed as a partial success. 

The presence or absence of a break in the specific star- 
formation rate as a function of stellar mass and the evolution 
of this with redshift appears to place significant constraints 
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Figure 11. Specific FIR luminosity as a function of redshift in stellar mass classes for "red" galaxies, i.e. those with SED types 1-2. 
Filled circles are the average of our estimates from the 70 &160u.m data which are the weighted average over four fields (Elais Nl, Elais 
N2, Lockman and CDFS). Thick vertical error bars come from the field-to-field scatter and thin error bars with hats from the variation 
between 70 & 160u.m estimates. Solid lines are our power-law trends with parameters give in Table[3] The dashed line shows the inverse 
Hubble time, galaxies above this line are 'bursting' producing stars at a higher rate than their historical average. 



on feedback models. Our large sample means we have small 
statistical errors and so present results in finer redshift and 
mass bins. This means our results are less subject to changes 
across a bin. Our work extends previous work to lower red- 
shift and has reduced errors at higher mass. By using the 
longer Spitzer wavelength and comparing 70 and 160n.m we 
provide a better constraint on the bolometric FIR power. 
Our large area means we cover a representative range of envi- 
ronments (so are not subject to sampling variance problems) 
and are better able to assess systematic errors by compar- 
ing independent fields. However, we have to compare with 
deeper samples to probe lower mass galaxies. It is worth 
emphasising that the stacking technique can overcome the 
problems of low sensitivity of far infrared data. Given the 
availability of wide area far infrared and submm surveys 
from Spitzer and in the future from Herschel and SCUBA-2 
a homogeneous analysis over representative samples of lower 
stellar mass would be possible with deeper optical data. 



6 CONCLUSION 

• We have used an improved stacking analysis to probe 
the far-infrared emission of galaxies near the peak of the 
cosmic infrared background. 



• We have shown that SWIRE sources with 5 , 3.6u.m con- 
tribute 70-80% of the cosmic infrared background with the 
main uncertainty in the determination of background itself. 

• We show that about 50% of the CIRB can be explored 
with these techniques to r < 23.5, arguing that deeper op- 
tical data will have a dramatic impact on the science that 
can be done in these fields. 

• We have measured the average specific FIR luminosity 
or specific star formation rate as a function of stellar mass 
and redshift. 

• We have found a trend sSFR oc M* with p ~ 
—0.38. This contrasts with SAMs in which the specific star- 
formation rate is constant with mass until a dramatic drop 
at high mass. This trend is stronger for early type galaxies 
(P ~ -0.46) than late type galaxies (/3 0.15). 

• We have found a strong evolutionary trend sSFR oc 
(1 + z) a with a = 4.4 ± 0.3 for 10.5 < log 10 M,/M Q < 12, 
steeper than a semi-analytic model which has a ~ 2.4. . 

• For early type galaxies the average evolution in this 
mass range is stronger (a ~ 5.7) but decreases to higher 
mass. 

• For late-type galaxies the trend is weaker but appar- 
ently independent of stellar mass, giving a mean rate a = 
3.36 ±0.16. 
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